Quantum size effects in Pb layers with absorbed Kondo adatoms: Determination of 

the exchange coupling constant 
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We consider the magnetic interaction of manganese phtalocyanine (MnPc) absorbed on Pb layers 
which were grown on a Si substrate. We perform an ab-initio calculation of the density of states and 
Kondo temperature as function of the number of Pb monolayers. Comparison to experimental data 
[Phys. Rev. Lett. 99, 256601 (2007)] then allows us to determine the exchange coupling constant J 
between the spins of the adsorbed molecules and those of the Pb host. This approach gives rise to 
a general and reliable method for obtaining J, by combining experimental and numerical results. 



Strong correlation phenomena built the fundament of 
a number of extraordinary effects observed in condensed 
matter systems. One of the most drastic manifestations 
is the Kondo effect^ arising from strong coupling between 
single localized spins and the conduction electrons. For 
an antiferromagnetic-type interaction below the charac- 
teristic Kondo temperature, Tk, this coupling cannot be 
treated perturbatively. It comes along with the forma- 
tion of a narrow peak in the density of states (DOS) of 
a Kondo adatom pinned to the Fermi energy of the con- 
duction electrons, £p. Its lineshape can be described by 
the well-known Doniach-Sunjic formula^ 3 -. 

The appearance of the Kondo resonance has various 
physical consequences. In the context of the transport 
properties of bulk diluted magnetic alloys it manifests as 
enhancement of the resistance at low temperatures. In 
mesoscopic systems consisting of a single quantum dot 
coupled to a pair of one-dimensional leads^ it results in 
a drastic enhancement of the conductance below Tr^JL. 

A system which recently has attracted much attention, 
both from the experimental and theoretical side, is com- 
posed of a metallic host surface with a Kondo adatom 
(adsorbed atom) and a STM tip, non- or ferromagnetic 9 . 
Contrasting the case of a single quantum dot coupled to 
leads, such a system demonstrates a much greater variety 
of transport regimes. This arises from the possibility of 
quantum interference between distinct transport paths, 
namely direct tip-host tunneling and indirect tip-adatom- 
host tunneling-^. In particular, the latter process gets 
extremely efficient for temperatures below Tk due to the 
appearance of a Kondo resonance in the adatom spec- 
tral function. The quantum interference mechanism leads 
to an asymmetry of the Kondo resonance. Instead of a 
Lorentzian lineshape it acquires a Fano-type lineshape, 
which is characteristic of interference between a discret 
level and a continuous spectrum. The appearence of this 
kind of a Fano-Kondo resonance has been found in vari- 
ous experiments with nonmagnetic tip a 10 ! 11 ' 12 ^ 3 and has 
been discussed in a number of theoretical work a 15 i 16 ' 17 . 

In a recent experimental work the quantum size effect 
on the Kondo temperature was investigated by Y.-S. Fu 
et alii*. The system under study consisted of individual 
manganese phtalocyanine (MnPc) molecules absorbed on 
ultrathm Pb films grown on Si(lll). Scanning tunneling 
spectroscopy shows asymmetric Fano-Kondo peaks in the 
differential conductance. One central implication of the 
work by Y.-S. Fu et alJ^ is that the variation of the film 
thickness, given by the number of Pb monolayers (MLs), 



leads to drastic changes in Tk, more specifically to pro- 
nounced oscillations with a period of 2 MLs. The lowest 
(23 K) and highest (419 K) values occured at thicknesses 
of 15 and 17 MLs, respectively. Because the Kondo tem- 
perature scales with the electronic DOS at the Fermi en- 
ergy, p = DOS(e_f), according to 



(1) 



where A is an amplitude and J is the exchange coupling 
constant, this observation has been attributed to oscilla- 
tions in the Pb DOS with the film thickness. Indeed, for 
every 2 MLs increase of the thickness one empty quanti- 
zation band moves down below the Fermi level and gets 
occupied^. 

In the following we present a state-of-the-art ab-initio 
investigation of the electronic structure of Pb films on a 
Si(lll) surface. Determining p = DOS(ef) as a function 
of the film thickness and using Eq.[T]we can perform a fit 
to the experimental data given in Fig. 2 of Ref. [3- Be- 
sides the amplitude A, the exchange coupling constant 
J is the only fitting parameter and therefore can be ob- 
tained without any further assumption. 

Our analysis is based on the augmented spherical wave 
approach, applying density functional theory within the 
local density approximation 19 . This method is particu- 
larly suitable for dealing with unit cells containing many 
atom a 20 i 21 ' 22 , which are needed to describe interfaces, 
since it utilizes a minimal atomic-like basis set. For the 
present system this basis set consists of Si 3s, 3p, 3d as 
well as Pb 6s, 6p, bd states, and is complemented by or- 
bitals of additional augmentation spheres. Brillouin zone 
integrations are performed using the linear tetrahedron 
method, where we checked the convergence of the calcu- 
lation with respect the fineness of the k-mesh by means of 
a growing number of up to 231 k-points in the irreducible 
wedge of the supercell Brillouin zone. 

Our Si(lll)/Pb(lll) supercells comprise a Si slab of 
about 25 A thickness, an attached Pb slab ranging from 1 
to 18 Pb MLs, and a vacuum slab on top. For generating 
the supercells we have started from the hexagonal repre- 
sentation of the fee lattice, with the (111) direction along 
the hexagonal Ch cx axis. For the Pb slab we use the ex- 
perimental lattice constant apb = 4.95 A. Moreover, for 
attaching the Si slab (asi = 5.43 A) to the Pb slab we 
have to artificially shrink the former and introduce a lat- 
tice strain. Although this lattice strain in the Si slab is 
substantial, no drawback on the Pb electronic structure 
is to be expected since there are no Si states in the vicin- 
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FIG. 1: Top: Total DOS at e F for a Pb(lll) film on Si(lll), 
depicted as a function of the number of Pb MLs. The values 
are normalized by the numbers of Pb atoms in the supercells. 
Bottom: Zoom for the 11-18 MLs range. 

ity of the Fermi energy. Our supercells extend 121.25 A 
along Chox and have as basal plane a parallelogram with 
a 60° interior angle and edges of length apb- 

We obtain for the DOS at the Fermi energy the results 
displayed in Fig. [T] as a funtion of the Pb film thickness. 
One clearly identifies the previously discussed oscillatory 
behavior, with a period of 2 MLs, up to a thickness of 9 
MLs. Deviations from this regular scheme are found in 
the intervals 9-12 and 14-16 MLs where the oszillations 
are suppressed^. Comparing Fig. [1] to the experimental 
datai^ indicates that our curve is shifted by 1 ML to the 
left hand side, see particularly the position of the 14-16 
MLs anomaly. This difference might either trace back to 
the experimental difficulty in counting the MLs or to the 
effects of the (planar) MnPc molecule on the local DOS, 
which is not included in our treatment. Since there is no 
ambiguity in the choice of the matching conditions at the 
Si-Pb interface, our approach elsewise should not lead to 
a 1 ML shift. 

Using Eq. Q]we next calculate from our DOS data the 
dependence of the Kondo temperature on the film thick- 
ness. To that aim, we use the second highest (150 K at 12 
MLs) and smallest (23 K at 15 MLs) of the experimen- 
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FIG. 2: Kondo temperature as a function of the thickness of 
the Pb(lll) film, given by the number of Pb MLs. 



tally observed Tk values to fit the free parameters A and 
J. The correctness of the highest measured temperature 
(419 K at 15 MLs) is quite dubious, see the discussion of 
the authors themselves, and thus not taken into consid- 
eration. At the Fermi energy, the ab-initio DOS amounts 
to 0.318 eV" 1 for 12 MLs and to 0.267 eV" 1 for 15 MLs. 
As a consequence, we obtain from Eq.[T]for the exchange 
coupling constant 

J « 0.32 ±0.04 eV. 

As compared to semiconductor bulk systems^, the value 
of J is 2-4 times smaller, which is not surprising because 
the coupling to a surface in general is much weaker than 
the coupling within the bulk. 

Due to the logarithmic dependence on Tk, the deter- 
mined value of J is particularly insensitive to experimen- 
tal inaccuracies. The outcome of our fitting procedure is 
displayed in Fig. [2l As to be expected, the dependence of 
Tk on the film thickness largely resembles the behaviour 
of p, see the bottom of Fig.[TJ We do not find indications 
of a dependence of the exchange coupling constant J on 
the film thickness between 11 and 18 MLs. Hence, there 
is a qualitative agreement with the experimental results, 
even though clarification of the 1 ML shift and Tk spike 
at 17 MLs should be tackled by future experiments. 

In summary, we have performed first-principles calcu- 
lations to establish the density of states of Pb(lll) films 
grown on a Si(lll) substrate. These ab-initio data have 
enabled us to fit experimental results from measurements 
addressing the dependence of the Kondo temperature on 
the Pb film thickness. In particular, we have succeeded 
in calculating the exchange coupling constant J, which is 
not accessible to direct experimental determination. For 
the system under investigation, a value of about J = 1/3 
eV transpires. 
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